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Abstract Fractals are a basic tool to phenomenologically describe natural objects having 
a high degree of temporal or spatial variability. From a physical point of view the fractal 
properties of natural systems can also be interpreted by using the standard formalism of 
thermodynamical fluctuations. In the present paper we introduce and analyze the fractal 
dimension of the Intra - Day- Variability (IDV) light-curves of the BL Lac objects, in 
the optical, radio and X-ray bands, respectively. A general description of the fluctuation 
spectrum of these systems based on general thermodynamical principles is also proposed. 
Based on the general fractal properties of a given physical system, we also introduce the 
predictability index for the IDV light curves. We have explicitly determined the fractal 
dimension for the R-band observations of five blazars, as well as for the radio band ob- 
servations of the compact extragalactic radio source J 1128H-5925, and of several X-ray 
sources. Our results show that the fractal dimension of the optical and X-ray observations 
indicates an almost pure "Brownian noise" (random walk) spectrum, with a very low pre- 
dictability index, while in the radio band the predictability index is much higher We have 
also studied the spectral properties of the IDV light curves, and we have shown that their 
spectral index is very closely correlated with the corresponding fractal dimension. 
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1 INTRODUCTION 

Fracta ls are a basic tool to phenomenologically describe natural objects. According to [Mandelbrot 
(Il983h . a fractal is a set whose Hausdorff dimension is not an integer The fractal dimension is one possi- 
ble parameter that characterizes chaotic systems, and the analysis of time series is one of the most com- 
mon means to find the fractal dimension from observable quantities. Fractals are also very useful for the 
analysis of the waveforms, a term that describes the shape of a wave, usually drawn as insta ntaneous val - 
ues of a variable quantity versus time. The fractal analysis of waveforms was introduced bv lKatzl(ll988h . 
who con sidered the p ossibility that the complexity of a waveform may be represented by its fractal di- 
mension. |Kat3 lUla) proposed that the fractal dimension can be measured empirically by sampling the 
waveform at N points evenly spaced on the abscissa, which discretizes the waveform into A^' = — 1 
segments. Then the fractal dimension can be obtained as 1?^? = \ogN' / [log A^' + log (d/L)], where 
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d = max [dist(i, j)] is the planar extent of the curve, and L = X^iLo dist(i, i + 1) is the length of the 
curve, where "max" stands for the maximum dist(i, j), the distance between points i and j of the curve. 
An alternative and mor e efficient algorithm for the computation of the fractal dimension was proposed 
bv LSevcikl(ll998U2006h . which we will consider in the following. 

The fractal analysis has been exte nsively used in astron omy and astrophysics for the study of the 
fractal properties of gamma ray bursts (IShakura et al.l 1994t . the critical properties of spherically sym- 
metric accretion in a fractal medium (iRoy feRaylbOOTlK, for th e analysis of fractal structures in th e 
photospheric and the coronal magnetic field of the Sun ( loshpa e t al. 2008; Di mitropoulou et ai]|2009l) , 
as a measure of the scale of homogeneity (lYadav et all 2010l) , for the study of the star forming re- 
gions (ISanchez et alj 2010l), an d for the analysis of the dark matter and gas distributions in the Mare- 
Nostrum universe ( Gaite 2010h , respectively. The quasar distribution on the celestial sphere is charac- 
terized by power laws as well wi th correlation dimension valu e equal from 1.49 to 1.58 for different 
redshift layers in the same range (iRozgacheva & Agapovll20TTI) . A fractal cosmological model, which 
ac counts for the observable fract al properties of the large-scale structure of the Universe, was discussed 
in lRozgacheva & Agapovl(l201ll) . 

Many classes of Active Galactic Nuclei have variability properties. The time scale of variability can 
be classified as yearly, monthly or daily. If the sources vary within one night or day, the time variation 
is called "Intra-day Variability", or simply IDV. There is one main method for confirming if the object 
is with IDV in the differential light curve, or not. This is the so-called statistical method. In this method 
one checks if the variability ampli tude is gre ater than 3 times the sigma value, or not. The best criteria is 
however 5 tim es of s igma criteria (" FaiilbOllI) . Another method, based on the so called C-parameter, was 
introduced bv lRomer o et al. (1999). The criteria for confirming the object is with IDV is the "C" value 
greater than or equal to 2.576. If this is the case, the IDV variability is confirmed at 99% confidence 
level. Furthermore, the object is said to be variable if the "C" value is gre ater than or equal to 2.5 76 at 
two different bands. The IDV can be checked by the formula proposed bv lHeidt & WagneJ (Il996l) . For 
checking the IDV for the source in radio band a statistical parameter called modu lation index is used . 
Usually, for the IDV light curve of blazars there is no significant pattern and period. iRoland et al.l (l2009b 
proposed that the formation and the rotation of a warp in the inner part of the accretion disk can lead to 
some perturbations of the beam that finally produce the observed IDV. 

It is the purpose of the present paper to analyze the IDV light curves by using their fractal dimen- 
sion properties, and to try to interpret the actual value of the fractal dimension by using the general 
formalism of thermodynamic fluctuations. The actual IDV signal is always fluctuating. Profiles of IDV 
signals are also remarkably varied. The possible mechanisms for generating IDV emissions could also 
be different. The fractal dimension could be provide a hint into the physical processes that could trig- 
ger the IDV processes. In the present pape r we have explicit ly determined the fractal dimension for 
the R-band optical emission of five blasars dGupta et al. 20081) and fo r the radio band emission of the 
compact extragalactic radio source J 1 128-H5925 (iGabanvi et al.ll2007l) . The fractal dimension aUows us 
the introduction of the predictability index, which represents a powerful indicator of the nature of the 
IDV signal. The general fractal properties of the IDV signals can be interpreted in the framework of the 
general theory of thermodynamic fluctuations for self-gravitating systems. The fractal dimension of the 
optical observations indicates an almost pure "Brownian noise" (random walk) spectrum, with a very 
low predictability index, while in the radio band the predictability index is much higher We have also 
studied the spectral properties of the IDV light curves by using the periodogram method, and we have 
shown that their spectral index is very closely correlated with the corresponding fractal dimension. 

The present paper is organized as follows. In Section|2]we consider the theory of the thermodynamic 
fluctuations of the self-gravitating systems, and obtain the general theoretical predictions for the fractal 
dimension of the fluctuation spectrum. The methods for the calculation of the fractal dimension and of 
the power spectrum are presented in Section |3] The fractal dimensions and the power spectra of several 
IDV signals in optical, radio and X-ray bands are obtained in Section H) The correlation between the 
fractal dimension and the spectral index is also discussed in detail. We discuss and conclude our results 
in Section|5] 
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2 FRACTAL PROPERTIES OF FLUCTUATING SELF-GRAVITATING SYSTEMS 

In the following we consider a system of self-gravitating particles (an astrophysical accretion disk, for 
example), embedded in a heat bath (an external medium) of gravitationally non-interacting particles. 
The system interacts with the surrounding bath through a friction force, due to which the particles in 
Brownian motion lose energy to the medium, and simultaneously gains energy from the random kicks 
of the thermal bath, which can be descr ibed by a stochastic force (for a description of the dynamics of 
the system via a Langevin equation see lLeung et al.l ( 1201 lb ). As a result of the interaction between the 
system and the heat bath, the thermodynamical quantities of the system fluctuate. In the following we 
consider the standard thermodynamic fluctuation formalism (iLandau & Lifshit3l 19801: iKandrupIl 19831) 
as applied to self-gravitating systems, and we derive the time dependence of the fluctuation spectrum. 

According to the fundamental principles of statistical mechanics, the probability w for a quantity 
to have a value in the interval from r to r + dr is proportional to e'^'^*'^'"'', where ASt{r) is the change 
in the entropy during the fluctuation (iLandau & Lifshit3ll980l) . The entropy change can be written as 
ASt = —Wi-nin/T, where VFmin is the minimum work needed to carry out reversibly the given change 
in the thermodynamics quantities in a small part of the system. Thus w cx e~^"""/^. For Wmin we can 
adopt the expression 

Wrr^ir^^AE-TAS + PAV + ARr^l^, (1) 

where AE, AS, AV are the fluctuations in energy, entropy, and volume, respectively, while Ai?inin is 
the minimal work necessary for reversible removal of mass AM for a distance Sr in the gravitational 
field of a mass M. Thus we have 

/ AE-TAS + PAV Ai?„,i„\ 
w oc exp ( — 1 . (2) 

Eq. ^ can be applied to any fluctuation, small or large. However, in the foUowing we consider only the 
case of small fluctuations. In this we can expand AE is series, obtaining 



AE=^AS+'-^AV+'- 
dS dV 2 



(ASf + 2^^ASAV + — (Avf 



(3) 



By taking into account that ^ = T, ^ = — P, we obtain 



AE ~TAS + PAV = i 

1 
2 



d^E _ „ d^E . d^E 



i^Sf + 2—-ASAV+ — iAVf 



dS^ ' ' dSdV dV^ 
{AS AT - APAV) , (4) 



givmg 



'ASAT-APAV Ai?,„, 
« cxp ( — — ] . (5) 



By taking V and T as independent variables, we obtain first 

AP=f^l AT+f^l AV. ,7) 



dTJy \dVyj, 



where C„ is the specific heat at constant volume. Therefore for the fluctuation probability we obtain 



w oc exp 



(8) 
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In order to es timate i?min we consider that it can be obtained as jRozgacheva & Agapovll2oTll 
lRozgachevalll984 

Ai?mi„ « G^^^Ar = ^GprAMAr, (9) 

3 

where M = Anpr^ /3 is the mass at radius r, and AM — Anpr'^Ar is the mass of a spherical shell 
of thi ckness Ar. Therefore the probability of thermal fl uctuations in a self-gravitating system takes the 
form (iRozgacheva & Agapovl20"TltlRozgachevall984 









1 


w cx exp 


24 




2(52 



Ar 

r 



where 5% ~ C., ^, and 



f 127rp 



1 



kT 



dp_ 

dV 



SttGp 



(10) 



(11) 



where k is Boltzman's constant. The root- sum-square relative density fluctuation (fluctuation spectrum) 
is given by dRozgacheva & Agapovl201ll) 




= 3 




dP 
dV 



^-i _l 



-1/2 



The fluctuations of the density are related to the fluctuations of the radial coordinate by 

Ap ^ _^Ar 



(12) 



(13) 



Since our main goal is to study the time variability of some astrophysical phenomena, we have to convert 
the radial variability into time variability. There are two possible time-scales for astrophysical phenom- 
ena, namely, the oscillation time scale and the free-fall time scale. For a homogeneous uniform self- 
gravitating system that behaves as a perfect gas throughout, the pulsation period tp is tp = 2t:/ ^J'inGp, 
while the free-fall time is tff — ^y4^T/3Gp. In both cases the time-scale is proportional to 1/ y/p. 
Therefore, assuming such a general time - density dependence, we obtain 

At 1 Ap 



With the use of Eq. ( fTsT i we obtain 



givmg 



t 

2 At 

sT 

r (X t 



2 P 

Ar 

r 

2/3 



(14) 

(15) 
(16) 



We can now discuss the limiting cases of the time-spectrum in a self-gravitating system. At spatial scales 
when the pressure gradients dominate, we have 




J Uwp 
\ kT 



dP 
dV 



-1/2 



oc r 



.-3/2 



cx t 



(17) 



At large scales, when gravitational effects are important, we obtain 

r 




-5/2 



oc t 



-5/3 



(18) 
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However, the above discussion is incomplete in the sense that in a typical fluctuating system we have 
two independent modes, one growing, and o ne decaying. A bet ter estimation of the time dependence 
of the fluctuations can be obtained as follows dRozgacheval 19881) . According to Chebyshev's inequality 
from the pro bability theory, the probability that a density fluctuation IS.pl p exceeds (5 > is given by 
(lFelledll97lh 

P(|Ap/p|><5)<^, (19) 

where C is a constant. Now let's assume that for the independent growing mode we have 8\ ^ t", 
while the decaying mode behaves as 82 ^ , where a, 7 > are constants. Then from Chebyshev's 
inequality we find Pi (|Ap/p| = 5i) ~ t~^", and P2 (jAp/pl = ^2) ~ t^^ , respectively. We now take 
into account the collective effects, and look for the probability for the existence of a growing density 
contrast if there are neighboring regions of space in which the density inhomogeneity decreases. From 
a mathematical point of view we need to find the conditional probability P12, which, according to the 
rules of probability theory, is given by Pi = Pi 2^2- Therefore 



P(|Ap/p| = <5) = Pi2 = ^~^-'"-'^. 



By using again Chebyshev's inequality we find 

Ap 



(20) 



(21) 



3 FRACTAL DIMENSION AND PREDICTABILITY INDEX 



The fractal dimension ca n be defined by the self-similar power law scaling function (IMandelbrotll 1 983l 
IRangarajan & Santll2004 

y{x)^asX^,D>0, (22) 

where y{x) denotes the number of self-similar objects in the sphere or circle of a radius x; Qs and D 
stand for the scaling factor and for the spatial fractal dimension, respectively. Physical systems evolve 
not only in space but also in time. Therefore many time dependent natural processes can be successfully 
fitted by the temporal counterpart of the fractal function given by Eq. (l22l i. 



y{t) = att°',t>0, 



(23) 



in which y{t) characterizes the time-evolution of the system, Dt is its temporal fractal dimension 
whereas at is a scaling factor. 

In the following we will obtain an expressio n for the cal culation of the fractal dimension by starting 
from the definition of the Haussdorf dimension (ISevciklll998>,2006.) . The Haussdorf Dh dimension of 
a set in a metric space is defined as 



D 



H 



lim 



ln[iV(e)] 
hie ' 



(24) 



where N (e) is the number of open balls of a radius e needed to cover the entire set. In a metric space, 
given any point P, an open ball of center P and radius e, is a set of all points x for which dist (P, x) < e. 
A line of length L may be divided into N{e) = L/ (2e) segments of length 2e, and may be covered by 
open balls of radius e. Thus, Eq. ( l24l i may be rewritten as 



Dh = lim 



- ln(P) + In (2e) 



In (6) 



lim 



1 - 



ln(L) - In (2) 



= lim 


[i 




[ In (6) J 



(25) 



Waveforms are planar curves in a space with coordinates usually having different units. Since the 
topology of a metric space does not change under linear transformation, it is convenient to linearly 
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transform a waveform into another in a normalized space, where all axes are equal. LSevcikI (Il998h 
proposed to use two linear transformations that map the original waveform into another, embedded in 
an equivalent metric space. The first transformation, normalizes every point in the abscissa as x* = 
Xi/xi-nax, where Xj are the original values of the abscissa, and .Tmax is the maximum Xj. The second 
transformation normalizes the ordinate as y* = {iji — y,nax) / (2/max ~ Vmin)^ where ?/,; are the original 
values of the ordinate, and ymin and ^max are the minimum and maximum y^, respectively. These two 
linear transformations map the N points of the waveform into another that belongs to a unit square. This 
unit square may be visualized as covered by a grid of N x N cells, N of them containing one point 
of the transformed waveform. Calculating L of the transformed waveform and taking e = 1/(2 x N'), 
where N' = N - 1, Eq. ^ becomes (Sevcik 199S) 

= Mi) ^ (26) 

" ln(2x7V') ' 

The approximation for the calculation of Dh improves as A^' — > oo. 

If the fractal dimension D for the time series is 1.5, there is no correlation between amplitude 
changes corresponding to two successive time intervals. Therefore, no trend in amplitude can be dis- 
cerned from the time series, and hence the process is unpredictable. However, as the fractal dimension 
decreases to 1, the process becomes more and more predictable as it exhibits "persistence". That is, 
the future trend is more and more likely to follow an established trend. As the fractal dimension in- 
creases from 1.5 to 2, the process exhibits "anti-persistence". That is, a decrease in the amplitude of 
the process is more likely to lead to an increase in the future. Hence, the predictability again increases 
dRangarajan & Sa nt 200 4,). 

Based on the fractal dimension we can introduce the predictability index P of a given physical 
system, defined as dRangarajan & San3ll997l) . 

P/ = 2|L»-1.5|, (27) 

where | | denotes the absolute value of the number D. We use absolute values since predictability 
increases in both the following cases - when the fractal dimension becomes less than 1.5 and when it 
becomes greater than 1 .5. In the former case, we have correlation (persistence) behavior and in the latter 
case, anti-correlation (anti-persistence) behavior. However, in either case, the process becomes more 
predictable. Thus, use of absolute values ensures that a process with D = 1.3 has the same predictability 
index as a process with D = 1.7. 



4 FRACTAL DIMENSION, PREDICTABILITY INDEX AND POWER SPECTRUM INDEX 
OF IDV LIGHT CURVES 

In order to analyze the fractal properties of IDV light curves, we consider a sample of several blazars, 
which ha ve been observed i n a simultaneous multi-wavelength observing campaign in their outburst 
phase bv iGupta et al.l (l2008b . Optical photometric monitoring of nine blazars was carried out in 13 
observing nights during the observing run of 2006 October 27-2007 March 20, using the 1 .02 m optical 
telescope equipped with CCD detector and BVRI Johnson broadband filters at Yunnan Astronomical 
Observatory, Kunming, China. The R-band observations of five blazars are represented in Figs. [T]and|2] 
respectively. 

There are many ways to characterize different noise sources . One possibility, which we have also 
explored in the present paper, is to consider the spectral density, that is, the mean square fluctuation 
at any particular frequency /, and its variation with frequency. In many natural processes the power 
spectrum P{f) is proportional to 1/ f^, where the power index /3 > 0. When /3 = 0, the noise is 
referred to as wh ite noise, when it is 2 it is referred to as a Brownian noise, and when it is 1 it is called 
as the 1/f noise dHavUn et al.lll995l:lDenisovlll998l) . 

In order to ob tain t he power spectr al index (3 we have used the periodogram method, as developed 
in IVaughanI (l2005b and lVaughanI (|20IO|) . respectively. To fit the data we have used power law fitting for 



Fractal properties of IDV light curves 



7 



all our bands. In fact, we can also use other options. However, since the number of available data points 
in the optical and radio bands is limited, for the sake of uniqueness we can only use power law fitting 
method. Although the X-ray data can be fitted by different functions, the optical data and the radio data 
cannot be fitted by other methods, such as the bending power method. Therefore, we fitted the data in all 
bands by using a power law. With the use of the above procedures, we can obtain the power spectrum, 
the power spectrum index, and all the other useful parameter information, such as the p- value. The error 
values after fitting have also been obtained, and we present them in Table [T] and |2l respectively, as a 
percentage error. 

The fractal dimensions, the predictability index PIj dv and the power index of the observed blazar 
optical time series are represented in Table [T] 



Object 


Date of observation 


D 


PIlDV 




Fit error of 13 


S2 0109+224 


2007-01-11 


1.52591 


0.051 


0.909647 


0.1723201 


S5 0716+714 


2007-01-11 


1.44996 


0.100 


1.41473 


0.1404088 


S5 0716+714 


2007-01-12 


1.45699 


0.086 


1.389167 


0.07811043 


S5 0716+714 


2007-02-23 


1.35088 


0.298 


1.672849 


0.2028525 


S5 0716+714 


2007-03-19 


1.57989 


0.159 


0.2951006 


0.1540101 


S5 0716+714 


2007-03-20 


1.44484 


0.110 


1.21574 


0.08695967 


PKS 0735+178 


2007-01-11 


1.56557 


0.131 


0.6411014 


0.1548084 


ON 231 


2007-01-11 


1.45285 


0.094 


0.7989009 


0.1536858 


lES 2344+514 


2007-01-12 


1.42142 


0.157 


1.148679 


0.2762177 



Table 1 Fra ctal dim ensions, predictability indices, and power indices of a sample of R-band 
light curves (iGupta et al...2008l) . 



As a second example of the determination of the fractal dimension we consider the short time-scale 
radio variations of the compact extragalactic radio source J 1 128+5925 jGabanvi et al. 2007). The flux 
density variability of J 1128+5925 was monitored with dense time sampling between 2.7 and 10.45 
GHz. For our analysis we use the data observed with the 100 m Effelsberg radio telescope of the Max- 
Planck-Institut fur Radioastronomie at 2.70 GHz, 4.85 GHz, and 10.45 GHz. Ten observing sessions 
that lasted several days during the period between 2004 - 2006 were performed. 

Apart from the optical and radio data, we also try to use short time-scale of IDV data from XMM 
database for calculating the fractal dimension. We reduced the data and followed the result from the 
published targets as shown in lGaur et al.l (l2010l) . From the XMM database we have reduced 7 blazars. 

We used XMM-Newton Science Analysis System (SAS) version 11.0 to reduce the EPIC-pn data 
taken from the archive. We followed the same procedure outUned in Gaur et al. (2010) to extract the 
lightcurves. In brief, we excluded all the high background due to flaring, and restricted all the events 
to the 0.3-10 keV energy band. We used 45 arcsec radius as the extraction region for both the source 
and the background (from source-free region). In addition, the background subtracted light curves were 
corrected for vignetting, bad pixels. Point Spread Function (PSF) variation, quantum efficiency, and 
dead time. All the light curves were binned with 100-second timing resolution. 

The radio and X-ray data, together with the corresponding power spectra are shown in Fig. |3] while 
the fractal dimension, predictability indices and the power spectral indices are presented in Table 2. 

The re is a relation between the fractal dimension D of the Brownian noise and the power index /3, 
namely (iHavUn et al.lll995l:lDenisovlll998l) 

D-^^■ (28) 

If the fractal dimension of the considered IDV sources is around D = 1.5, then the power spectral index 
of the spectra is around /3 = 2. We have studied the correlation between the fractal dimension D and 
the power spectral index of the optical and X-ray IDV sources - we did not include in our study the 
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Object 


Date ot observation 


D 


PllDV 


1^ 


Fit error ot j3 


J 1128+5925 


2005-09-16 


1.28642 


0.42716 


2.188192 


0.1472077 


J 1128+5925 


2006-04-28 


1.25587 


0.48826 


2.042647 


0.1100852 


ON 231 


2002-06-26 


1.43877 


0.122 


1.621563 


0.06751574 


lES 2344+514 


2005-06-19 


1.60232 


0.204 


0.63486 


0.06634579 


PKS2155-304 


2002-05-24 


1.44012 


0.11976 


1.055890 


0.0554563 


PKS2155-304 


2002-05-24 


1.56946 


0.13892 


0.6418146 


0.06857538 


PKS2155-304 


2002-05-24 


1.33378 


0.33244 


1.615948 


0.06764645 


PKS2155-304 


2000-05-31 


1.57255 


0.1451 


1.023354 


0.04338111 


PKS2155-304 


2000-11-21 


1.51078 


0.U2156 


1.171058 


0.04761672 



Table 2 Fractal dimensions, predictabilit y indices, and power spectral indices of the compact 
extragalactic radio source J 1128+5925 dOabanvi et ailbOO?!) an d of the X -ray sources ON 



231, lES 2344+514 and PKS 2155-304 respectively (|GauretalJ|2010|) 



radio data, due to their small number. There is a very good correlation between D and /3, and for this 15 
sources the relation between the fractal dimension and the power index can be given as 



D = 1.66229- 0.16699 x /?. 
Eq. ( |29] l can also be written, with a very good approximation, as 

10-/3 



D 
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(29) 



(30) 



The correlation statistics of the D — /3 correlation is represented in Table |3] For the chosen set of 
data the correlation coefficient between these two parameters is 0.8428. 



c. c. 


C.C.^ 


Standard Error 


Total number of cases 


0.8428 


0.7104 


0.0469 


15 



Table 3 Regression statistics for the correlation between fractal dimensions and power spec- 
tral indices of optical and X-ray sources. The correlation coefficient is denoted by C.C. 



The statistical significance of the correlation coefficient is presented in Table |4] by means of the 
ANOVA analysis. In Table[3] C.C represents the correlation coefficient. It represents the fitting accuracy. 
With the C.C equal to 1 and -1 the scatter diagram can be fitted by one single straight line. C.C^ 
represents the square of C.C. In Statistics, the number of degrees of freedom is the number of values 
in the final calculation of a statistic that are free to vary. If the F-statistics computed in the ANOVA 
table is less than the F-table statistics, or the p-value is greater than the alpha level of significance, then 
there is no reason to reject the null hypothesis that all the means are the same. For our p-level value, 
we obtained 0.00008, and thus we conclude that the null hypothesis should not be rejected at the 1% 
significance level, since the p-level value is much less than 0.01. 

The relation between the theoretical data fit, given by Eq. (|29] | and the observational data is pre- 
sented graphically in Fig.|4] The presence of the strong Brownian noise in the IDV spectra also strongly 
support the ide a that random stochastic processes play a significant role in their formation and evolution, 
as suggested in iLeung et al.l (l20Ilh . 

5 DISCUSSIONS AND FINAL REMARKS 

Since the determination of the fractal dimension of the times series is relatively simple, it can provide 
useful information about the signal fluctuation from variable astrophysical sources. In the present paper 
we have explicitly obtained the fractal dimensions for several IDV emissions in both optical and radio 
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Name 


d.t. 




MS 


F 


p-level 


Regression 


1 


0.07041 


0.07041 


31.89023 


0.00008 


Residual 


13 


0.0287 


0.00221 






Total 


14 


0.09912 









Table 4 ANOVA Table for the correlation coefficient of the fractal dimensions and power 
spectral indices of optical and X-ray sources, d.f. represents the number of degrees of free- 
dom, SS is the sum of squares, MS is the mean sum squares value, F represents the value of 
the F-statistics, and p-level represents the measure of the degree of the insignificance of the 
regression model. 



bands. From the point of view of the fractal analysis the optical and the radio signals shows some re- 
markable differences. In the optical and X-ray bands the fractal dimension of the IDV signals is very 
close to the value D = 1.5, indicating an almost perfect "Brownian noise" (random walk, o r Wiener pro- 
cess) fluctuation spec trum (Brownian noise can be generated by integrating white noise) jHavlin et al.l 
ll995l : lDenisov|[l998l) . The predictability index of this signal is very low, showing that the emission in 
the optical band is dominated by purely stochastic processes. On the other hand the radio band emission 
pattern of the compact extragalactic radio source J 1128H-5925 shows a relatively low fractal dimen- 
sion of the order of D = 1.25, with a much higher predictability index. This result can be interpreted 
from a physical point of view in the framework of the thermodynamic theory of fluctuations as follows. 
Assuming that in both fluctuating phases the self-gravitating fluctuations have a time dependence given 
by Eq. ( fT6] l. we have a = 7 = 2/3. Then from Eq. (ISTT l it follows that the time dependence of the fluc- 
tuation spectrum is of the order of ~ t°'+t ^ ^4/3 ^ ^1-33^ which gives a fractal dimension D = 1.33, 
very close to the observed values of D ^ 1.25 and D = 1.28, respectively. Therefore this implies that 
the dominant time scale in the emission process is of the order of the gravitational free fall time scale. 
On the other hand, in order to obtain through the same method the fractal dimension {D = 1.5) in 
the optical phase of the IDV the assumption of a t^O-'i's dependence of the fluctuation spectrum in the 
fluctuating phases is necessary. This shows that the characteristic time scale of the density fluctuations 
cannot be of the same order of magnitude as the free fall time, and other physical processes, like, for 
example, viscous dissipation, may be responsible for the dynamics. 

We have also obtained the general time dependence of the fluctuation spectrum in a self-gravitating 
system, as given by Eqs. ( [TtI i and ( fTsT i, respectively. The density fluctuations can be related to the 
luminosity fluctuations of the emitting source. By assuming that the emission of the radiation is in the 
form of thermal radiation, from the Stefan-Boltzmann law we find that the luminosity (energy emitted 
per second) from a spherical surface with radius R is given hy L = AnE^aT'^j where a is Stefan's 
constant, and Te// is the effective temperature of the emitting surface. With the use of the hydrostatic 
equilibrium equation and of the virial theorem one can show that L oc (IChoudhurill2"010l) . With 
the use of this relation we obtain immediately /S.L/L cx ^pj p- Therefore, with the use of Eqs. (fTTI i 
and ( fTSl l the radial distance and time spectrum of the fluctuations of the luminosity in self-gravitating 
systems can be written as 



\\\Lj 



ocr-^^' ocr'^K (31) 



The constants Dr and Dt can be interpreted physically as the fractal dimensions, describing the 
fractal properties of the thermal radiation emission spectra from fluctuating sources. The comparison of 
the fractal dimension of these spectra with the IDV observational data could give a powerful indication 
about the nature of the physical processes (thermodynamical or gravitational, random or deterministic) 
that play the dominant role in the electromagnetic emission process. 

In the present paper we have also studied the correlation between the fractal dimension and the 
power spectral index /3. We have found a very good correlation between these two physical parameters. 
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Fractal properties, or self-similarity, can be found in many aspects in nature. Our analysis has con- 
centrated on the fractal properties of the IDV signal. Due to the limited amount of data used, our results 
can be considered only as a first step in the investigation of the fractal properties of the IDV's. It will be 
important in the future to investigate more signals at different wavelengths, in order to obtain a more ac- 
curate determination of the fractal dimension, and to clarify the relation between the fractal dimensions 
of different bands and the spectral power index, respectively. 
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APPENDIX 

To obtain the spectral index /3 and for fitting our data, we used the periodogram method as developed 
in Vaughan (2010). The periodogram is used for dealing wit h the evenly sa mpled time series Xk of K 
points at intervals AT. The definition of the periodogram is (IVaughanl2005l) 

2AT 

where P is the modulus-squared of the discrete Fourier transform. 

The unit of the normalization of the periodogram is (rms/mcan)^ Hz^^. The initial purpose of the 
periodogram was to search for "hidden periodicity" in time series analysis. But, for a single time series, 
the periodogram of a noise process shows a great deal of scatter around the underlying power spectrum. 
For the case of a given frequency, / (fj) is scattered around the true power spectrum, P (fj). It also 
follows a distribution with two degrees of freedom, 

Hfj)^P{fj)xl/2 (33) 

where is a random variable distributed as with two degrees of freedom. It is equivalent with an 
exponential probability density, with a mean and variance of two and four, respectively, 

p^2{x)=e^p{^x/2)/2. (34) 

For the case of white noise flat spectrum, the power level is calculated by making use of Eqs. ( [33] ) 
and ( l34b to obtain the maximum likelihood of a given periodogram. Similarly, for the more general case 
of non-white noise, we need to examine the spurious peaks in the periodogram, and regard it as Brownian 
(red) noise. A simple method to obtain a reasonable estimate of power-law spectrum is S{f) — af~^. 
a is the so called normalization of power-law power spectrum, and /3 is the power spectral index. 
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S5 0716+714, 2007-02-23 




Fig. 1 R-band variability and power spectram of the sample of blazars observed in lGupta et al.l 
(l2QQ8h . 
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Fig. 2 R-band variability and power spectrum of the sample of blazars observed in lGupta et al.l 
(l2()08l) . 
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Fig. 3 Radio- band variability and power spectram of the compact extragalactic radio source 
J 1 128+5925 (IGabanyi et al.ll2()07h and X-ray variability and power-spectrum of the sources 
ON 231 and lES 2344+514, respectively. 
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